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Spacelike and timelike response of confined relativistic particles 
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Abstract. Basic theoretical issues relating to the response of confined relativistic particles are considered 
including the scaling of the response in spacelike and timelike regions of momentum transfer and the role 
of final state interactions. A simple single particle potential model incorporating relativity and linear con- 
finement is solved exactly and its response is calculated. The response is studied in common approximation 
schemes and it is found that final state interactions efi'ects persist in the limit that the three-momentum 
transferred to the target is large. The fact that the particles are bound leads to a non-zero response in the 
timelike region of four-momentum transfer equal to about 10% of the total strength. The strength in the 
timelike region must be taken into account to fulfill the particle number sum rule. 
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1 introduction 

Deep inelastic scattering (DIS) of leptons by hadrons is 
generally discussed in the framework of the naive parton 
model and the QCD-improved parton model using the 
operator product expansion. This approach has been 
very successful in determining the evolution of the struc- 
ture functions as a function of the square of the four- 
momentum transferred to the hadron.Q In the leading 
order of the model the hadron is approximated by a col- 
lection of noninteracting quarks and giuons. The struck 
quark is assumed to be on the mass-shell both before and 
after its interaction with the electron. Basic theoretical 
considerations bring the validity of these assumptions into 
question. Q 

Based on the assumption that the struck constituent 
is on the mass-shell before and after interaction with the 
probe, the response is predicted to be in the spacelike 
region for which the energy transfer i' is less than the 
magnitude of momentum transfer, |q|, as a consequence 
of the inequality. 



i^=J\k + q\' 



— 



(1) 



In the many-body theory (MBT) one expects, at least 
naively, that final state interactions (FSI) should have 
an effect on inclusive scattering cross sections with elec- 
tromagnetic probes from systems whose constituents are 
confined. Scattering of high energy probes from compos- 
ite systems, such as electron scattering by nuclei [Band 
nucleons , or neutron scattering by liquid helium [|| , is 
often used to study the structure of the bound system. 
The common assumption is that in DIS at sufficiently 
high energy the probe is incoherently scattered by the 
constituents of the system. In the plane wave impulse ap- 
proximation (PWIA), which neglects FSI effects, DIS is 
directly related to the momentum and energy distribution 
of the constituents in the target. 

The role of FSI effects has been studied extensively 
in electron scattering from nuclear targets [^,|^ and neu- 
tron scattering from liquid helium |^ . Recently it has been 
suggested that they may also influence DIS of leptons by 
hadrons In the present study we focus on scattering 
from targets with confined constituents. The correspond- 
ing physical case concerns DIS from nucleons where, in dis- 
tinction from the nuclear and liquid helium cases, the con- 
stituents are confined in both the initial and final states. 



Here k and are the momentum and mass of the struck 
quark, respectively. The predicted response is discontinu- 
ous at the boundary |q| = between space and timelike 
regions. In fact interactions among the constituents in the 
initial state take the constituents off the mass-shell and 
move response of the target into the timelike region of 
four-momentum transfer. 
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2 The response and scaling variables 

We consider the response to a hypothetical scalar probe 
which couples to the density of a single scalar constituent. 
This allows us to ignore complications due to spin and the 
Lorentz structure of the response though it retains the 
qualitative features of a more realistic model where one 
considers the coupling of a spin-i fermion to the conserved 
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Fig. 1. The jql-i' plane. The spacelike region is above the 
|q| = V line and timelike is below. Lines of constant Q'^ > 
are parabolas which lie entirely in the spacelike region and 
approach = v a& v ^ oo. The observed (g^ > 0) response 
of the proton lies in the shaded area. 



electromagnetic current. The response is 
i?(q,z.)=^|(/|^e''i--^|0)p5(ii;,- 

/ .7 



Eo - ly) (2) 



where J2j is over all the particles and the J^i over all 
energy eigenstates. It is viewed as the distribution of the 
strength of the state s*'* ""^ |0) over the energy eigen- 
states of the system having momentum q. It is not neces- 
sarily zero in the timelike, > |q| region. 



2.1 Scaling variables 

The conventional variables of the parton model, = 
|qp — and the Bjorken x — Q"^ /2Mi>, used to describe 
the DIS structure functions of a hadron of mass M, are 
confined to the spacelike region of the |q|-J^ plane for pos- 
itive values of accessible in lepton scattering experi- 
ments, as shown in Fig. |. The observed (Q^ > q) DIS 
response is limited to a narrow region in the |q| — plane 
illustrated in Fig. ^ It is bounded by the elastic limit, 
Uei = i/lqp + M^ — M on one side, and by the photon 
line on the other. In the limit of large |q| the width of the 
observed response at fixed |q| is M . Lines of constant 
intersect the elastic limit curve at x = 1 and approach the 
photon line at small x. 

We wish to study the full range of response possible 
for a system of bound constituents including the region 
of timelike momentum transfer. Therefore we study the 



response, R{(\,u) as a function of v and |q| in the rest 
frame of the system as is common practice in the MBT. 
Lines of constant |q| in Fig. cross the photon line [v — 
|q|) and go into the timelike region. The natural scaling 
variable in the MBT approach to DIS |^ is y = i^ — |q|. At 
large |q| the response is expected to depend only on y, and 
not on q and v independently. This variable is equivalent 
to the Nachtmann variable ^ since 



e = ^(|q| 



(3) 



In the limit of large the — x, thus y scaling includes 
Bjorken scaling. However, both y and ^ span both space- 
like and timelike regions at fixed |q| unlike x at fixed . 



3 Model calculation 

We have studied the exact response of a simple "toy" 
model which contains the basic features of relativity and 
confinement to obtain further insights on the possible re- 
sponse in the timelike region and it's effects on the sum 
rules. In this model we assume that the response of the 
hadron is due to a single light valence quark confined 
within the hadron by its interaction with an infinitely mas- 
sive color charge. We model this interaction by a linear 
fiux-tube potential, and use the single particle Hamilto- 
nian, 

H = y^|p|2 + m2 + Vct r (4) 

containing the relativistic kinetic energy operator. In the 
limit mq = used here, the H can be cast in the form: 



V|2 



(5) 



where p' = p/cr^/'*, and r' = cr^/^r are dimensionless. The 
response i?(|q|, v) of the model then depends only on the 
dimensionless variables |q'| — jql/cr^/^ and v' = v/a^^^. 
The main conclusions of this work are independent of the 
assumed value of ct; however, we show results in familiar 
units using the typical value ^/a = 1 GeV/fm. 

The model may be viewed as that of a meson with a 
heavy antiquark or that of a baryon with a heavy diquark. 
It is obviously too simple to address the observed response 
of hadrons. For example, it omits the sea quarks and ra- 
diative gluon effects contained in the DGLAP equations 
to describe scaling violations. Nevertheless its exact 
solutions are interesting and useful to study scaling, the 
approach to scaling, and the contribution of the timelike 
region to sum rules. A similar model has been considered 
by Isgur et al. 

The Hamiltonian is diagonalized in the spherical mo- 
mentum basis and the response is calculated to ensure 
that the full strength of the integrated response. 



i?(q, iy)di' ~ 1, 



(6) 



is obtained in the chosen basis for all values of the mo- 
mentum transfer considered in this work with < 0.02 % 
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error. In order to obtain a smooth response we assume 
decay widths for all the excited states dependent on the 
excitation energy v. 

Figure ^ shows the response calculated for values of 
|q| > 3 GeV as a function of y. The scaling behavior 
is clearly exhibited; at large |q| the i?(|q|,j/) becomes a 
function /(y) alone. This scaling is equivalent to ^ scaling 
via (Eq. ||). 

In Fig. H we show the response at various values of 
|q| < 2 GeV compared with that for |q| = 10 GeV, to 
study the approach to scaling. At small |q| the scattering 
is dominated by resonances, and the first inelastic peak is 
due to the lowest excited state with n = 1 and £ — I, 335 
MeV above the ground state. In our toy model, the elastic 
scattering occurs at = or y = — |q|, since our hadron is 
heavy. This elastic scattering contribution is omitted from 
Fig. I ^ 

For y ^ 0, i.e. for small ^, the response approximately 
scales at relatively small values of |q|, comparable to cr^/^. 
As |q| increases, the range over which scaling occurs is ex- 
tended to more negative values of y, i.e. to larger values of 
^. The contribution of each resonance shifts to lower y and 
decreases in magnitude following the -R(|q| — > oo,y). This 
behavior is seen in the experimental data on the proton 
and deuteron [ p^ and interpreted as evidence for quark- 
hadron duality. Thus the toy model seems to describe 
some of the observed properties of the DIS response of 
nucleons. It exhibits y or equivalently ^ scaling at large 
|q| as observed and an approach to ^ scaling similar 
to that seen in recent experiments. 



= ^(0|e'i-('"^-'"')|0). (7) 

When q is large only the i — j terms in the above sum 
contribute, and therefore the integral gives the number 
of particles in the system. In contrast the sums of the 
response in the parton model are obtained by integrating 
the response over ^ > at fixed Q^. These sums will 
fulfill the particle number sum rule only if the response 
in the timelike region is zero. As mentioned earlier, the 
response of a collection of noninteracting particles lies in 
the spacelike region. Interaction effects, however, can shift 
a part of the strength to the timelike region. Evidence for 
shifts caused by interactions is discussed in Ref. 

Returning to the "toy" model the i?(|q|, and there- 
fore the f{y) extend into the timelike {y > 0) region. The 
sum- rule given by Eq. (||) , counts the number of particles 
in the target. It is necessary to integrate over the timelike 
region to fulfill this sum rule. About 10% of the sum is 
in that region independent of ^/a. The response expressed 
as R{Q^,S^) also scales at large where |q| is necessar- 
ily large. It becomes a function of ^ alone. However, the 
integral: 

poo p\<i\ 

/ i?(Q2 ^ oo, £,)d^ = / i?(|q| oo, iy)diy < 0.9, 
Jq Jo 

(8) 

because the contribution of the timelike region is omitted. 
Here we have defined ^ = |q| — without the conventional 
1/M scale [Eq.(|)]. 



3.1 Particle number sum rule 

In general the particle number sum rule in MET is ob- 
tained by integrating the response at large |q| over all 
v>0: 

poo 



4 Final state interaction effects 

We study the effects of the FSI of the struck particle on 
the response. Analytic calculations of the width of the 
response are presented for a general spherically symmet- 
ric potential and numerical results for a linear confin- 
ing potential are given. These indicate that the FSI in- 
crease the width of the response beyond that predicted by 
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PWIA. The analytic calculations also consider the nonrel- 
ativistic problem, in which q is large compared to aU the 
momenta in the target, but smaller than the constituent 
mass m. The main differences between the nonrelativistic 
and the relativistic response are that the former peaks at 
V = |qp/2m and has a width proportional to |q|, while 
the latter peaks at v ^ |q|, and has a constant width in 
the scaling limit. 



4.1 Moments of the response 

In the case of a single confined particle, the state of the 
system after the probe has struck the target is 



\X) 



|0), 



(9) 



where |0) denotes the ground state of the particle. The 
state \X) is not an eigenstate of the Hamiltonian and 
therefore has a distribution in energy. It has a unit norm, 



{X\X) 



|g-»q rg»q r|Q<) ^ ^ rpj^g ^^^-g^j strength of the 



response, given by the static structure function 



Si\q\)= / d^^i?(|qU), 
Jo 



(10) 



is therefore unity. In many-body systems 5'(|q|) is not nec- 
essarily equal to one. Subsequent formulas pertain to the 
general case and show factors of 5'(|q|) cxphcitly. 

The mean excitation energy of the state \X) is given 
by the first moment of the response: 



1 



1 



^(|q|) 



sm Jo 



{X\H~Eo\X) = ^j^ I di^i^Ri\q\,iy). 

(11) 

The width of the distribution in energy is characterized 
by the second moment of the energy about the mean: 



Z\2( 



1 , ,/ (X\H\X)\\ , 



(12) 



Substitution of the Eq. (||) into the formulas for the first 
three moments of the response give the following results: 

^^(|q|) = |q| + (V')o - i^O 

A'i\ci\)^^{k^)o + (V\~(V)l 

+ ^((fc'^>o-(fc^Wo)+o(^) (14) 

Here O (^j^^ denotes the neglected terms of that and 

higher order and the angle brackets with subscript '0' in- 
dicate averaging with respect to the ground state. Thus 
'^(IqI) = IqI ~ {T)o in the limit |q| oo, where T = 
IpP is the kinetic energy. The requirement that F(|q|) — 
|q| becomes constant is naturally satisfied in this limit. 
These expression demonstrates that the average energy 



and width of the exact response is independent of |q| in the 
limit |q| — > oo as necessary for y scaling. It also shows that 
the width has a kinematic contribution dependent upon 
the target momentum distribution, and an additional in- 
teraction contribution. 

As mentioned, the PWIA assumes that a constituent 
of momentum k, after being struck by the probe, may be 
described by a plane wave with momentum k -f q in an 
assumed average potential chosen to give the exact 17 of 
Eq.(n3h. From the PWIA response we calculate 



^PWIA 



1 



o 



1 



(15) 



contains only the first term of the exact result [Eq.([l4|)] 
due to the target momentum distribution. The second 
term, {V^)o - {V)l of Eq.(|l|) represents the FSI contri- 
bution neglected in the PWIA. It does not vanish in the 
|q| ^ oo limit for relativistic kinematics. 

|2 

In the non-relativistic case, i/jvj?, = + Vif), the 
exact T> is given by: 



vnr 



J-l 

2to' 



(16) 



For the width of the NR-PWIA response we obtain: 



^'NR.-PWIAi\^\) 



1 



+ ~ {kYo) ■ (17) 

Note that in Eqs.(|6l) and (|l|) we have not taken the 
|q| oo limit. 

The width of the exact NR response is: 



^^fl(|q|; 



— ^NR~PWIA 

+ {V^)o - {V)l 



+ -{{Ve)o-{V)o{k^)o) 



(18) 



It differs from Anb-pwia in terms of order l/|q| which 
can be neglected in the scaling limit. Thus, in contrast to 
the relativistic case, the FSI do not increase the width of 
the NR-PWIA response at large |q|. 

Finally we consider the on-shell approximation (OSA) 
in which the energy of the struck constituent is that of 
a free relativistic particle before and after the interaction 
with probe, as assumed in the quark-parton model. The 
response in OSA depends only on the momentum distri- 
bution of target constituents and obeys y scaling. The 
average excitation in OSA is 



77o5^ = |q|-(T)o + ^(fc2)o + o(^ 
and the width is given by: 



(19) 



(20) 



The exact value of v [Eq.(|l3|)] is reproduced by the OSA 
for any potential. However, the A'qqj^ has (fc^)o — {T)'o in 
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Fig. 4. The response versus y calculated exactly for Jo = 100 
MeV (thin solid curve) and Fo = 50 MeV (dotted curve). The 
response in OSA are shown for |q| = 10 GeV (dashed) and 
|qt oo (dot-dashed). The PWIA response for |q| = 10 GeV 
and |q| oo lie on essentially the same (thick solid) curve. 



place of the {V'^)q — (V)'^ in the leading term of the exact 
[Eq. (|l^] . For a massless particle in a linear confining 
potential, i.e. for the Hamiltonian of Eq.(||), (T)o = {V)o, 
and (fc^)o = {V'^)o. Therefore for this particular Hamilto- 
nian the OSA reproduces the exact value of Z\; but the 
shape is wrong. 



4.2 Numerical results 

We first compare the response functions for |q| = 10 GeV 
before comparing their moments. In Ref. ||lj| it has been 
shown that the scaling limit is obtained for such values of 
|q|. The exact response, Eq.(||), is a sequence of 5 functions 
at v = Ei — Eq. In order to obtain a smooth response we 
assume decay widths Iq for all the excited states. Note 
that the energies of the states |/) that contribute to the 
response at |q| = 10 GeV are large, therefore their decay 
widths are not affected by the energy dependent terms 
assumed in Ref. 1 14 . The response including decay widths 
is given by: 



The OSA results in the discontinuous curves shown in 
Fig. ^. They are discontinuous at the lightline (|q| = v) 
because the response of free particles is limited to the 
spacelike region ly < |q|. The discontinuity at y = is in 
clear conflict with the exact response which is continuous 
across the lightline and is non-zero in the timelike (y > 0) 
region. Therefore the OSA appears to be unsatisfactory 
even though for the special case of a linear potential it 
has the exact values of S'(|q|), i^(|q|) and Z\(|q|). 



References 

1. R. K. Ellis, W. J. Stirling, and B. R. Webber, QCD and 
Collider Physics (Cambridge University Press, Cambridge, 
1996), pg. 108. 

2. G. Altarelli and G. Parisi, Nuc. Phys. B 126, 298 (1977). 

3. J. D. Bjorken, in 7th Conference on the Intersections of 
Particle and Nuclear P hysics (AIP Pre ss, Quebec City, 
Quebec, Canada, 2000), [hep-th/0008044 

4. B. Frois and I. Sick, Modem topics in electron scattering 
(World Scientific, Singapore, 1991). 

5. R. N. Silver and P. E. Sokol, Momentum Distributions 
(Plenum Press, New York, 1989). 

6. O. Benhar et ai, Phys. Rev. C 44, 2328 (1991). 

7. O. Benhar and V. R. Pandharipande, Phys. Rev. C 47, 
2218 (1993). 

8. S. J. Brodsky et ai, Phys. Rev. D 65, 114025 (2002). 

9. O. Benhar, V. R. Pandharipande, and I. Sick, Phys. Lett. 
B489, 131 (2000). 

10. O. Nachtmann, Nucl. Phys. B 63, 237 (1973). 

11. R. L. Jafi'e, Lectures presented at the 1985 Los Alamos 
School on Quark Nuclear Physics, 1985. 

12. N. Isgur, S. Jeschonnek, W. Melnitchouk, and J. W. Van 
Orden, Phys. Rev. D 64, 054005 (2001). 

13. L Niculescu et ai, Phys. Rev. Lett. 85, 1182 (2000). 

14. M. W. Paris and V. R. Pandharipande, Phys. Lett. B514, 
361 (2001). 



^(q,i^) = y"l(^|e''' '"|o)P ( ^ . 

^ Y V27ry {Ei-Eo-i^y+rg/i 

(21) 

The responses obtained with Fq ~ 100 and 50 MeV are 
shown in Fig. |4| along with the PWIA and OSA responses 
for |q| = 10 GeV and for |q| — > oo. The difference between 
the exact responses for Fq — 100 and 50 MeV are much 
smaller than those between the exact and the approxi- 
mate. 

We note that the shape of the PWIA response is qual- 
itatively similar to that of the exact, however, its width 
is too small. This is a direct consequence of the neglect of 
interaction terms in A [Eq.(p^] as discussed in the last 
section. The width A of the response is 409 MeV, while 
the ApwiA = 326 MeV. 



